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In this work we reports results of ferromagnetic resonance studies of a 6% 15 nm (Ga,Mn) As 
layer, deposited on (OOl)-oriented GaAs. The measurements were performed with in-plane oriented 
magnetic field, in the temperature range between 5 K and 120 K. We observe a temperature induced 
reorientation of the effective in- plane easy axis from [TlO] to [110] direction close to the Curie 
temperature. The behavior of magnetization is described by anisotropy fields, H e g (= 4nM — H2±), 
H 2 \\ , and . In order to precisely investigate this reorientation, numerical values of anisotropy fields 
have been determined using powerful - but still largely unknown - interval calculations. In simulation 
mode this approach makes possible to find all the resonance fields for arbitrarily oriented sample, 
which is generally intractable analytically. In 'fitting' mode we effectively utilize full experimental 
information, not only those measurements performed in special, distinguished directions, to reliably 
estimate the values of important physical parameters as well as their uncertainties and correlations. 

PACS numbers: 07.05.Kf, 68.47.Fg, 75.30. Gw, 75.50.Pp, 75.70.-i, 75.70.Ak, 75.70.Cn, 75.70.Rf, 76.50.+g 



I. MOTIVATION 

Despite numerous and intensive studies, an origin of 
the in-plane uniaxial magnetic anisotropy in (Ga,Mn)As 
remains unknown. However, both its strength and orien- 
tation can be described on the ground of p-d Zener model 
assuming the existence of a fictitious epitaxial strain 
PHI]. On the other hand, as this is the leading mag- 
netic anisotropy at elevated temperatures and that the 
means of its control have already been demonstrated [1] , 
further studies on this intriguing property are timely and 
important. In particular, a presence of the temperature- 
induced 90° rotation of the direction of the easy axis 
[l[ may ease the vector manipulation of magnetization 
in future devices. In this communication we report on 
the technical analysis and results of FMR studies of such 
a (Ga,Mn)As layer that exhibit the easy axis rotation at 
temperatures close to its Curie temperature. 

The free energy density for our system, expressed by 
anisotropy fields (-ff's) rather than by more customary 
anisotropy constants (-f^'s), has the form (gusHM (Zee- 
man), and inactive out-of-plane fourfold anisotropy term 
have been omitted): 

F = 2 7T M 2 cos 2 9 (1) 
MH 2± cos 2 9 (2) 
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where ([T]) is shape anisotropy, © - ordinary uniaxial out- 
of-plane anisotropy, ([3]) - uniaxial in-plane anisotropy, 
and ([4]) - fourfold, in-plane anisotropy. 



Here polar angles angles 9 and tp refer to the orienta- 
tion of the magnetization vector, M, not to the orienta- 
tion of an external field H . 

Original experimental data, taken at fixed frequency 
9.378 GHz, are shown in Fig. [TJ together with simulated 
spectra. 
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FIG. 1. Experimental values of FMR fields vs. external field 
orientation and temperature (points). Lines show computed 
spectra. 

The numerical values of three parameters: H c g = 
AttM - H 2 ±, (where H 2 ±_ = 2K 2J _/M), H 2 \\ and ff 4 || 
were determined for each temperature separately. 



II. OUTLINE OF NUMERICAL PROCEDURE 

We are using interval calculations not only for accurate 
simulations of FMR spectra when the values of all rele- 
vant parameters are known, but also to estimate ('fit') 
such parameters, together with their uncertainties and 
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correlations, utilizing full information hidden in experi- 
mental data. An excellent introduction to interval arith- 
metics can be found in [5|. Here we limit ourselves to 
the very brief description: i) an interval x is a bounded 
set of real numbers: x = [a, b] = {R 9 i : a ^ i O}, 

ii) it is possible to perform arithmetic operations, eval- 
uate functions, etc., using intervals instead of numbers, 

iii) the so obtained results are intervals as well, iv) the 
multidimensional intervals are called boxes for obvious 
reasons, and v) ordinary real numbers may be identified 
with ('thin') intervals, for example 7 = [7,7]. Similarly, 
/(x) is an interval containing all the possible results of 
evaluation of f(x) when x G x. Unfortunately, /(x) usu- 
ally overestimates the range of true values of f(x) - but 
always contains them all. 

The unconventionality of our approach to fixed fre- 
quency FMR data fitting is that we try to adjust un- 
known parameters in such a way that the classical for- 
mula 
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is satisfied for each experimental datum. 

The main difficulty lies in the fact that the partial 
derivatives of free energy density have to be evaluated 
at (stable) equilibrium position of magnetization vector, 
characterized by two (initially unknown) polar angles tp 
and 8, and being the solution(s) of the equation's system 



VF = (dF/dO, dF/dtp) = 0. 



(6) 



This makes the inversion of formula ([3]), to obtain res- 
onance field(s) as a function of microwave frequency w, 
essentially impossible. 

Our algorithm operates on a list C of boxes. At the 
beginning, the list contains only one member - the initial 
search domain. Further we perform following steps: 

1. select the biggest box from the list C 

2. bisect its longest edge obtaining two offspring 
boxes, then remove the parent box from the list C, 

3. investigate each one of the two offspring and: 

- discard it, if infeasible, or 

- put it back on the list C, otherwise. 

The above procedure is repeated until the list contains 
only small boxes. The rest is easy, provided the final clus- 
ter of boxes consists of only one connected component. 
For details on how to calculate mean values, variances 
and correlations between the searched parameters see @ . 
Here we only clarify when the box is considered infeasi- 
ble. First of all, each box is characterized by it's quality 
factor, computed as an interval quantity Q — [Q, Q]. For 
each resonance field (and its corresponding orientation) 
we try to evaluate the interval value of r.h.s. of expres- 
sion ([5]). If the upper bound of so calculated interval is 



negative, then the currently considered box of searched 
parameters is infeasible as it cannot describe any reso- 
nance at all. If all experimental data pass this test, then 
the ranges of their resonance frequencies can be found 
from ([5]) and the interval Q can be evaluated. Again, for 
each data element i the absolute difference between both 
sides of (J5J) is calculated, producing the interval Aw;. Fi- 
nally, Q = maxi (other choices are possible but we 
prefer this one). The box is considered infeasible when 
it's Q exceeds the reference value, which is equal to the 

lowest Q ever seen during calculations. 

It remains to comment on resonance frequency cal- 
culation. The angles {9, <p) are unknown and should be 
determined prior to resonance frequency evaluation, for 
each data element separately. Starting from full ignorance 
(8 G [0,7r], ip G [0,27r]), we divide this initial 2D box un- 
til its edges are shorter than, say, 0.5°. Discarded are all 
boxes satisfying either of two conditions: ^ dF/88 or 
^ dF/dip - they certainly cannot contain the equilib- 
rium position of the vector M. Needles to say that fail- 
ure to find such a position immediately invalidates the 
searched parameter box. 



III. RESULTS AND DISCUSSION 

Following the procedure described in previous section, 
we have computed the values of parameters determin- 
ing the free energy density (HHJ). During computation 
g was kept fixed at 2.00. The results are presented in 
Fig. [2] The numerical values are consistent with those 
reported by others 0-E1 , obtained by the same technique 
but at lower temperatures. Using those results, with- 
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FIG. 2. Anisotropy fields vs. temperature - computed val- 
ues. Lines are eye-guides only. Apparent lack of smoothness 
may be attributed to inaccurate temperature readings (the 
measurements were taken during two separate sessions). For 
better visibility, the curves Hi\\ and Hm are shifted horizon- 
tally by — 0.2K and +0.2 K, respectively. 

out any prior smoothing, we were able to determine the 
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directions of spontaneous magnetization in the interest- 
ing temperature range, near Tq . In Fig. [3] one can see 
that, in absence of the external field, there are only two 
such directions, antiparallel to each other, not four as one 
might expect from symmetry arguments. The presence 
of external magnetic field, oriented along [110] (Earth's 
field is sufficient) breaks even this symmetry, and the 
spontaneous magnetization aligns itself exactly along the 
external field above the transition temperature. Below 
the transition temperature the spontaneous magnetiza- 
tion deviates considerably from its 'natural' [110] (or, 
equivalently, [110]) position. This means that precise ex- 
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FIG. 3. Equilibrium position of magnetization vector vs. 
temperature for various strengths of the external field. The 
field is always directed along high-temperature easy axis [110], 
marked also as 135°, and its strength is as indicated. 



all should go to as T —> Tq - but they don't. In addi- 
tion, the simulated angular FMR dependencies (Fig. [T]) 
seem not so accurate as one might expect. Those intrigu- 
ing facts, together with high reliability of interval analy- 
sis, strongly suggest that the uniaxial in-plane anisotropy 
may be incorrectly accounted for in ([US]). Since the sam- 
ples exhibiting this special behavior are very thin, in nm 
range, then it is quite possible that the source of uniax- 
ial anisotropy may be related to surface effects. Indeed, 
AFM surface studies of MBE-grown samples [lOj revealed 
the presence of well ordered ripples, parallel to [110] di- 
rection. If this was true, then the free energy expression 
should be appended with the surface term [111 ], propor- 
tional to | cos (9? — 7r/4)|. This, however, is beyond the 
scope of a current paper. 

IV. CONCLUSIONS 

The interval calculus has been demonstrated to be 
a very powerful tool for difficult problems of experimen- 
tal data analysis. In particular, it is probably the only 
method able to utilize complete experimental informa- 
tion acquired during FMR measurements. Its unrivaled 
reliability allows us to state a sound hypothesis con- 
cerning the nature of the somewhat mysterious in-plane 
uniaxial magnetic anisotropy observed in thin layers of 
(Ga,Mn)As grown on (OOl)-oriented GaAs substrate. 

We have also shown that FMR technique is very help- 
ful in accurate tracing the temperature dependency of 
magnetic anisotropy in close vicinity of the Curie tem- 
perature, that is where the magnetometric data are least 
reliable. 



amination of spontaneous magnetization, especially of its 
components perpendicular to the external field, is diffi- 
cult. The temperature range, in which the reorientation 
occurs, is also very sensitive to the presence of even very 
weak field. 

The behavior of anisotropy fields near Trj, estimated 
as ~ 113 K for our sample, is somewhat unusual: they 
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